DNA breathing dynamics: Analytic results for distribution functions of relevant 
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We investigate DNA breathing dynamics by suggesting and examining several different Brow- 
nian functionals associated with bubble lifetime and reactivity. Bubble dynamics is described as 
an overdamped random walk in the number of broken base pairs. The walk takes place on the 
Poland-Scheraga free energy landscape. We suggest several probability distribution functions that 
characterize the breathing process, and adopt the recently studied backward Fokker-Planck method 
and the path decomposition method as elegant and flexible tools for deriving these distributions. In 
particular, for a bubble of an initial size xq, we derive analytical expressions for (i) the distribution 
P{tf\xo) of the first-passage time tf, characterizing the bubble lifetime, (ii) the distribution P{A\xo) 
of the area A till the first-passage time, providing information about the effective reactivity of the 
bubble to processes within the DNA, (iii) the distribution P{M) of the maximum bubble size M 
attained before the first-passage time, and (iv) the joint probability distribution P{M,tm) of the 
maximum bubble size M and the time tm of its occurrence before the first-passage time. These dis- 
tributions are analyzed in the limit of small and large bubble sizes. We supplement our analytical 
predictions with direct numerical simulations of the related Langevin equation, and obtain a very 
good agreement in the appropriate limits. The nontrivial scaling behavior of the various quantities 
analyzed here can, in principle, be explored experimentally. 

PACS numbers: 87.14.gk, 87.10.Mn, 02.50.-r, 05.40.-a 
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I. INTRODUCTION 

The Watson-Crick double helix structure of DNA de- 
rives its stability from the phosphodiester bonds in the 
single-stranded sugar backbone, and from the hydrogen 
bonds between complementary base pairs on opposite 
strands [H, Ql- In practice, access to the inside of the 
double helix, and therefore, the unzipping of a specific 
region of base pairs is essential for all physiological pro- 
cesses involving DNA, e.g., for replication, transcription, 
and protein binding |3| . 

Several mechanisms, like heating changing the pH 
of the environment and application of external force 
can lead to unzipping of the double-stranded DNA. 
This phenomenon is referred to as DNA denaturation. 
The process occurs progressively, starting with the dou- 
ble strand separating locally into single strands to form 
loops, or, "bubbles". These bubbles fluctuate in size 
through stepwise zipping and unzipping of the base pairs 
at the two zipper forks where the bubble connects to 
the double strand. At low temperatures, bubbles once 
formed eventually close again in time. With the increase 
of temperature, however, the bubbles grow in size in time 
to ultimately coalesce with neighboring bubbles and com- 
plete the denaturation process. The melting temperature 
Tm is defined as that at which half of the DNA molecule 
is denatured, and has typical values ^ 70 — 100°C for 
standard salt solutions [5[ . The number of bubbles varies 
from only a few ones well below T,„ up to several hun- 
dreds close to Tm- 

Breathing dynamics, referring to the dynamics of fluc- 



tuating DNA bubbles, has been a topic of intense re- 
search for many years 0, 0] ■ It has recently regained in- 
terest with the development of new experimental tools 
that allow for the direct observation of the dynamics 
of a single DNA molecule 0, On the theoretical 
side, various methods have been used to study differ- 
ent aspects of the breathing process d, [l^l , and to in- 
vestigate the interaction of the DNA with binding pro- 
teins: the master equation approach [TTI. [l^. a stochas- 
tic Gillespie scheme [l^ , the Fokker-Planck equation ap- 
proach based on the Poland-Scheraga free energy func- 
tion [l3 - [l6j , and stochastic dynamic simulations based 
on the Dauxios-Peyrard-Bishop model [13, Ell ■ Specifi- 
cally, the thermally-induced denaturation problem has 
been recently studied by mapping it onto a quantum 
Coulomb problem [l^, [20] . These studies have enhanced 
our understanding of general aspects of both polymer 
dynamics as well as specific biochemical processes. For 
recent advances, see [3- 

DNA breathing occurs on a timescale shorter than 
the equilibration time of the single strands forming the 
bubbles [1|. Based on this observation, breathing dy- 
namics may be regarded as a random walk in the one- 
dimensional coordinate x, the number of broken base 
pairs. Ignoring heterogeneity in the DNA structure, this 
random walk may be modeled as a noisy overdamped mo- 
tion at a finite temperature T on the Poland-Scheraga 
free energy landscape, J^{x) ^ jx + cksTlnx [l^, [Tst - 
The parameters c > and 7 (which can be of either 
sign) are defined later in the paper. As we show be- 
low, this form of the free energy implies a crossover scale 
Xch- For small bubbles {x < Xch), the random walk takes 
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place in a potential Inx. In the opposite limit, the 
potential grows linearly with x, thereby implying a dif- 
ferent dynamics. Furthermore, the sign of 7, as given by 
the system temperature, determines the nature (attrac- 
tive/repulsive) of the potential. For 7 > 0, which hap- 
pens at temperatures T < T™, the potential is attractive 
for all bubble sizes, thereby implying an eventual bubble 
closure (x ~ 0). On the other hand, above r,„, when 
7 < 0, a large bubble {x > Xch) evolves under a repulsive 
linear potential to grow in size toward full denaturation, 
while a small bubble (x < Xct) may still shrink in size to 
closure under the influence of the attractive In x poten- 
tial. 

In this work, following the above picture, we comple- 
ment previous single-bubble studies by suggesting and 
analyzing new measures for exploring the DNA breath- 
ing process. We focus on several first-passage "Brow- 
nian" functionals [2l| of the fluctuating bubble, which 
eventually closes again, and derive their probability dis- 
tribution functions (pdfs) . We separately study the small 
and large bubble limits, which exhibit different behav- 
iors. The functionals of interest characterize the lifetime 
of the bubble, the time-integrated bubble size till the 
first-passage time (the bubble "area"), its maximum size 
before closure, and the characteristic time for attaining 
the maximum size. These measures are relevant for es- 
timating the effective reactivity of the bubble, e.g., its 
efficiency for binding processes. 

Another objective of this work is to advocate the use 
of the recently studied backward Fokker-Planck (BFP) 
method [l^] and the path decomposition (PD) method 
[23| . which builds on the Feynman-Kac formalism p^ . 
for exploring DNA bubble dynamics. These techniques 
have been extremely useful in studying many aspects of 
classical Brownian motion, as well as for exploring related 
problems in computer science and astronomy ,22,„ 25., ^6] . 
Here, for the first time, we adopt these elegant methods 
in the context of DNA breathing dynamics. Using the 
BFP method, we derive and solve differential equations 
for the Laplace transforms of various Brownian function- 
als. This is in contrast to the standard Fokker-Planck 
treatment, which yields the distribution function to ob- 
tain a bubble of a given size at a given time [l^ - [l6l [isj . 
Utilizing the PD approach, we can calculate the distri- 
bution functions of interest by splitting a representative 
path of the dynamics into parts, and then considering 
the weight of each part separately. This is justified by 
the Markovian property of the dynamics. 

In order to gain a qualitative understanding of the 
DNA breathing process, we separately consider the cases 
of small and large bubbles. Our analysis reveals new 
scaling laws for the pdfs of various Brownian functionals, 
which are evidently distinct for small and large bubbles. 
We further compare our analytical predictions with di- 
rect numerical simulations of the corresponding Langevin 
equation and observe a very good agreement in the ap- 
propriate limits. 

The paper is organized as follows. In Section |lll we 



recall the random walk model and discuss the distribu- 
tion functions of interest and their relevance to the DNA 
dynamics. The BFP and the PD methods are also ex- 
plained, along with a short description of the numerical 
technique adopted. In Section lllll we study the dynamics 
of small bubbles, and derive the probability distribution 
functions of several first-passage Brownian functionals. 
In Section IIVI the dynamics of large bubbles is exam- 
ined. We draw our conclusions in Section FVl 



II. MODEL, QUANTITIES OF INTEREST, AND 
METHODS 

A. Model 

We follow the Poland-Scheraga approach, and inter- 
pret bubbles as occurring due to free-energy changes to 
the double-helical ground state Q. Measuring the size 
of a bubble by the number of broken base pairs, and de- 
noting this number by the continuous variable x > 0, the 
Poland-Scheraga free energy is given by Q 



J-{x) = 7o -t- 7a; -I- cksT In x, 



(1) 



where 70 is the free energy barrier to form the initial bub- 
ble, while the term 7a; stands for the free energy required 
in breaking x base pairs. The entropy loss in forming a 
closed polymer loop is taken into account by the term 
c/cbT In a;, where fc^ is the Boltzmann constant, T is the 
temperature, while c is a universal constant determined 
by the loop configurations |27|. In Eq. ([J), a cutoff at 



X ~ 1 is implied. The parameter 7 is assumed to have 
the simple form, 7 = 7i(T,„ — r)/T,„, where 71 — AksTr, 
with Tr = 310K being the reference temperature. 

At finite temperatures, the stochastic dynamics of 
DNA breathing can be modeled by the overdamped 
Langevin equation [ill , 



dx ^dF , 
dr dx 



(2) 



Here, ^(r) is a Gaussian white noise with (^(r)) = 0, and 
(C(t)C(t')) = 2DkBTS{T-T'). The kinetic coefficient D 
has the dimension of (A:bT)^^s^^. Using the free energy 
^ in Eq. and redefining the time variable, t = 

2DkBTT, results in the equation 



at X 
where Ci = c/2, C2 - 7i(r - T^) / {2kBTTm), and 

m) = 0, 

Equation ([T]) implies a crossover scale, 



(3) 



Ci 

\C2\ 



(4) 



(5) 
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such that for small bubbles, x < Xch, the free energy 
is dommated by the entropic term ^ lux. Correspond- 
ingly, the Langcvin dynamics ([3]) is essentially governed 
by the term —Ci/x. For large bubbles with x > Xch, 
the base-pair dissociation term ~ ^x dominates the free 
energy, and correspondingly, it is the C2 term which dic- 
tates the Langevin dynamics. For T < T^, the Langevin 
dynamics occurs in an attractive potential for all bubble 
sizes, thereby ensuring eventual bubble closure. Above 
Tjn, large bubbles with x > Xch grow in size under a 
repulsive linear potential to ultimately yield full denat- 
uration, while a small bubble with x < Xch niay evolve 
towards closure under the influence of the attractive In x 
potential. We will utilize below the length scale ([5]) in dis- 
tinguishing between small and large bubbles. Note that 
at precisely the melting temperature T™, when C2 = 
and the characteristic bubble size diverges, the Langevin 
dynamics becomes identical to that for small bubbles at 
all temperatures T ^ T^. 

Figure [T] depicts several characteristic paths evolving 
under Eq. ^ by adopting different values of the param- 
eter C2 with a fixed Ci. If IC2I > Ci and C2 < (i.e., 
T < Tm), the bubble closes sufficiently fast in compari- 
son to the case C2 = (top panel). In the opposite limit, 
taking positive values for C2, one observes the melting 
process reflected in the divergence of the first-passage 
time (bottom panel). 
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FIG. 1: (Color online) An illustration of several representative 
paths x{t) following the time evolution of Eq. ((3|. All the 
paths begin at xo = 4. The values of the parameters Ci and 
C2 are marked in the figure. These paths have been generated 
by using a Brownian simulation (see Section Til E|) . 



B. Quantities of interest 

Our primary focus is on several first-passage Brownian 
functionals of experimental relevance. We will consider 
the following quantities and explore their pdfs for small 
and large bubbles. 

(i) Bubble lifetime. The first-passage time pdf P{tf\xo) 



i.e., the pdf of the time of closure for bubbles of initial 
size Xq, provides information about bubble lifetime. A 
related quantity is the survival probability C{xo, t) = 1 — 
P{xo\tf)dtf which can be inferred from experiments 
by measuring fluorescence correlations of a tagged DNA 

m. 

(ii) Area under a path. For the path described by Eq. 
(|3]), we define the area under the path before the first- 
passage time as A = J^^ x{t')dt' , see Fig. [21 and calculate 
its pdf P{A\xo). This quantity is of interest since it pro- 
vides a measure for the effectiveness of breathing-assisted 
processes, i.e., binding of proteins to the reactive sites of 
the DNA bases. As an example, consider a process that 
can take place only inside the double helix, on the single- 
stranded DNA. Let us assume that this process is facili- 
tated with increasing bubble size, and that it requires a 
sufficiently long bubble lifetime. While the first-passage 
time distribution provides information about the average 
bubble lifetime, it does not contain any hint of the av- 
erage bubble size before closure. Similarly, P{M), the 
distribution of the maximum bubble opening before clo- 
sure provides a measure for the bubble size, but it does 
not inform us about the corresponding timescale. Thus, 
we propose the pdf P{A\xo) of the area A covered till 
the first-passage time as a useful quantity that provides 
a measure of bubble reactivity by containing information 
about both size of the bubble and its characteristic life- 
time. 

(iii) Maximum bubble size AI . Another proposed mea- 
sure for quantifying bubble reactivity is the distribution 
of the maximum bubble size before the first-passage time, 
P{M). Consider again a binding process taking place 
only inside the double helix. Assume next that, due to 
geometrical constraints, the process may materialize only 
when the bubble is large enough. If the timescale of this 
process is very short, shorter than the average bubble 
lifetime, a relevant measure for the bubble reactivity is 
its maximum opening before closure. 

Quantities (i) , (ii) and (iii) will be calculated below by 
following the backward Fokker-Planck method discussed 
in Section HTCl 

(iv) Maximum size M and the corresponding time tm- 
The joint probability distribution function P{M, tm) will 
be investigated here by following the PD method, which 
builds on the Feynman-Kac formalism [22|, [2^ ; see Sec- 
tion IIIDI Using this pdf, one can further calculate the 
distribution function P{tm) of the time at which the bub- 
ble attains its maximum size before closure. This latter 
pdf is of interest since it provides information about the 
(average) time of occurrence of the biggest bubble before 
closure. Processes taking place inside the DNA, facil- 
itated by increased bubble size, will most likely occur 
around that time. 

Figure [2] illustrates a typical path following Eq. ([3|) . 
The path begins at Xq and ends at the origin (bubble clo- 
sure), staying positive in between. The various measures 
suggested above arc also indicated in the figure. 
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FIG. 2: (Color online) An illustrative path x{t) that begins 
at xo and evolves under Eq. ((3|. Here, tf marks the time 
the path crosses the origin for the first time (corresponding to 
bubble closure), A is the area enclosed under the path, M and 
tm represent respectively the maximum value that the path 
reaches before the first-passage time, and the corresponding 
time of occurrence. This path has been generated by using a 
Brownian simulation (see Section fll E[) with Ci = 1, C2 = — 1, 
and xq = 2. 



C. The backward Fokker-Planck (BFP) method 

Following , we recall here how to calculate the sta- 
tistical properties of a Brownian functional, defined as 

T= [ ' U{x{T))dT. (6) 
Jo 

Here, x(t) is a path representing the motion ([3]) which 
starts at at time r = and propagates up to t = iy, 
the first-passage time. In the above equation, U{x{t)) 
is a specified function of the path, whose choice depends 
on the quantity of interest. For example, to compute 
the distribution of the first-passage time tf, one chooses 
U{x{t)) = 1. For the area distribution till the first- 
passage time, one should consider U{x{t)) = x. To find 
the pdf P(T|xo), noting that the random variable T can 
be only positive for these choices of ?7(a;(T)), one consid- 
ers its Laplace transform 

Jo 

where the angular brackets denote averaging over all 
paths starting at xq at r = and ending at the first 
time they cross the origin. For simplicity of notation, 
in what follows, we suppress the variable p in the func- 
tion Q(xo,p). In order to derive a differential equation 
for Q{xo), we follow [l^ and split the interval [0,^/] 
into two parts. During the first interval [0,Ar], the 
path starts from xq and propagates up to xq + Ax. In 



the second interval [At, i/], the path starts at Xq + Ax 
and reaches at tf. Here, At is a fixed, infinitesi- 
mally small time interval. We get, to leading order in 
At, J^^' U{x{T))dT w U{xo)AT + f^^ U{x)dT, and hence, 
from Eq. d?]), 

Q(x-o) « e-P^("«)^^(Q(a-o + Ax))A. 

w (l~pU(xn)AT){Q{x^+ Ax))ax. (8) 

Now, the average denoted by the angular brackets is 
performed over all realizations of Ax. The dynami- 
cal equation ([3]) gives Ax = i^(xo)AT -|- ^(0)Ar, with 
F(xo) = C2 - Ci/xq. Substituting for Ax in Eq. (P, 
expanding Q{xq + Ax) in powers of At, and averaging 
over the noise by using (^(0)) = and (^^(0)) = 1/At 
for small At, one obtains, to lowest order in At, the 
ordinary differential equation, 

l<mx^ ( C,\ dQjxo) 

(9) 

Boundary conditions. The above equation is valid for 
Xq G [0, 00] with the following boundary conditions: (i) 
For an infinitesimally small bubble, xo — > 0, the first 
passage time vanishes, tf — ?> 0, so that Q{xo = 0) = 1. 
(ii) If the bubble is initially large, xq — > 00. the first 
passage time diverges, hence, Q{xo — ?> 00) = 0. 

We emphasize that the differential equation re- 
ferred to as the backward Fokker-Planck equation [2^ . 
directly provides us with the Laplace-transformed pdfs 
of various quantities which are determined by the choice 
of U{x). In contrast, the standard Fokker-Planck method 
adopted in [l3 - [l6l [l8| yields the density distribution 
function P{x, t) to obtain a loop of size x at time t. Thus, 
these two approaches are distinct, providing complemen- 
tary information. 

D. The path decomposition (PD) method 

The principle of this technique is simple: Since the 
motion in Eq. ^ is Markovian, a typical path can be 
split into, e.g., two parts. Then, the weight of the whole 
|)a^h is the product of the weights of the two split parts 

The above idea allows us to calculate the joint proba- 
bility distribution P{M, tm) of the maximum bubble size 
M and the time tm at which this maximum occurs before 
closure, given that the initial size of the bubble is fixed 
at Xq G [0,A/]. By integrating over M, one can further 
obtain the marginal distribution P{tm)- Wc compute 
P{M,tm) by splitting a typical path into two parts, be- 
fore and after tm, with the respective weights Wl and 
Wr, so that the weight W of the whole path is 

W = WlX Wb.- (10) 

On the left side of tm, the path propagates from Xq aX t = 
to M — e at t = tm, without ever attaining the value 
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or M during the interval [0,tm] [111- The weight Wl can IIIIDI and IIVDI We explain here in some detail the cal- 



be determined by using a path integral treatment based 
on the Feynman-Kac formalism, as we explain below. On 
the right side of t^, the path starts from M — e at t = 
and ends at the origin at (with tf > tm), without 
crossing either the level M or the level in between. At 
the end of the calculation, one needs to take the limit 
0. 

The calculation of Wu will be explained in Sections 



dilation of ■ Since the white noise in Eq. ([3]) is Gaus- 
sian, the probability of a path is given by 



P[{x(r)}] cx exp 



1 



*■ /dx Ci ^ ^2 

— + C2 

\ar X 



(11) 



The weight Wl is then given as a sum over contributions 
from all possible paths. 



a;(t,„)=M- 



Wl oc 



'Dx{t) exp 



dr 



dx Ci N 2- 

-, 1 L'2 

cLt X 



X{e[x{T)]\{e[M-x{T)] 



T = 



T = 



Xq 



-M-e 



Ci 



a;(t,„)=J\/-e 



'Dx{t) exp 



(12) 



In the above equations, the terms rit=o ^'^'^ 
nt=o G[M — x{t)] enforce the requirements that the path 
does not cross either the level or the level M for times 
between and tm- Following Feynman-Kac [l^], the path 
integral in Eq. ([T2|) is identified with the propagator 
(M — eje"^*"* |a;o), corresponding to the quantum Hamil- 
tonian of a single particle of unit mass, 

+ (13) 
with h = 1. The potential energy V{x) is given by 

V{x) = \h{^-C-)' ifO<x<7\/, 

[00 if .T = or a; = M. 

Note that the infinite potential energy at a: = and at 
X = M enforces the requirement that the path never 
crosses either the level or the level M . Finally, we get 



E. Simulations 

The statistical properties of Brownian functionals 
studied here can be numerically obtained by integrating 
the overdamped Langcvin equation Using a second- 
order stochastic Runge-Kutta algorithm [2^, we update 
the trajectory by following the rule, 

x{M) = XQ + ]^[F{xa) + F{xa + F{xQ)M + To)]M 
+ro, (16) 

where F{x) = C2 — Ci/x. Here, Fq is a random number 
sampled from a Gaussian distribution with zero mean 
and width given by (Fg) = At. For all simulations pre- 
sented in this work, we take At = 10^"^, unless stated 
otherwise. We generate a large set of paths, all starting 
at a particular xq and ending close to the origin (within 
a preassigned numerical tolerance value) . Averaging over 
an ensemble, we generate various pdfs which we compare 
with our analytical results. 



c °° 
p=i 

(15) 

where ippix) and Ep are the eigenfunctions and eigenen- 
ergies, respectively, of the Hamiltonian H in Eq. (fT5|) . 
As a result of the infinite potential barrier at a; = and 
at a; = M, the eigenfunctions satisfy ipp{x = Q,M) — 0. 



III. SMALL BUBBLE DYNAMICS 

We begin our analysis by considering small bubbles, 
X < Xch, at all temperatures, T ^ Tm- The analysis is 
also valid for bubbles of all sizes at precisely the melting 
temperature Tm- In these cases, the nonlinear entropic 
term in the free energy ^ dictates the dynamics, result- 
ing in the Langevin equation. 
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where iiit)) = and = Ht - *')• The pdfs 

P{tf\xo) and P{A\xo) are obtained by the BFP method, 
where the differential equation that needs to be solved is 
given by Eq. © with C2 = 0, 



2" 



dxQ 



Ci dQjxo) 
xo dxo 



-pU{xo)Q{xo)=0. (18) 



The boundary conditions are (i) Q{xq — > 00) = 0, and 
(ii) Q{xo — > 0) = 1. We also derive analytical results for 
P{M) and P{M,tm), as explained below. 



A. First- passage time distribution: P{tf\xa) 

We compute the distribution oi tf, the time at which 
the bubble closes for the first time, assuming its initial 
size is fixed at xq, by substituting U{xo) — 1 in Eq. (ITSl) . 

Id^Qixo) 
2' 



Ci dQixo) 



dxl 



Xo dxo 

The general solution of Eq. (HH) is [s^ 



pQixo) = 0. (19) 



g(xo) = A/„(v^xo) + BK, 



2pxo 



(20) 



Here, Ia{x) and Ka{x) are the modified Bessel func- 
tions of the first and second kind, respectively. Also, 
a = Ci + 1/2, and A and B are arbitrary constants to be 
determined from the boundary conditions. Since for large 
X, Ia{x) ~ I\J2'KX and K^ix) ^ ^/^TJ2x [sij, in 
order to satisfy the condition Q{xo — > 00) ~ 0, we must 
have A — 0. To satisfy the condition Q{xo 0) = 1, 
we note that as x ^ 0, K^ix) « r(a)2"-Va;" for a > 
(all, which gives S = (y2^)'^i+i/V[r(Ci + 1/2)2^1-1/2]. 
Following these considerations, we get the particular so- 
lution. 



Q{xo) 



Ci+1/2 



(V2^) 



Ci+1/2 



r(Ci + 1/2)2^1 



On taking inverse Laplace transform, we get 



(21) 



^2Ci + l 
''0 



C,-3/2^-.l/2tf^ (22) 



r(Ci-f 1/2)2^1+1/2'^/ 

as obtained earlier in [l^ [s^ . 

Next, we compare the analytical prediction ([22]) with 
numerical simulations under the full bubble potential, 
including the C2 contribution, in order to explore the 
regime of validity of the above result, see the top panel 
of Fig. [31 The sampled trajectories all begin at Xq = 2. 
When Xch = C1/IC2I ^ 10 > xq; C2 < 0, we observe a 
good agreement between numerics and analytical results. 
For Xo ~ Xch, deviations occur since then the contribu- 
tion of C2 cannot be neglected (inset). 

Besides the distribution ([2^. other related quantities 
of experimental relevance are the moments, (i^), ob- 
tained from Eq. as 



xf r(Ci - fc + 1/2) 

2k r(Ci + 1/2) 



for fc > Ci + 1/2, (23) 



while (ty) diverges for k < Ci + 1/2. Another important 
quantity is the persistence, or survival probability of the 
bubble, defined as 



(24) 



C(xo,t) = l- / P{tf\xo)dtf 



where P{tf\xo)dtf sums up the probabilities of all 
events where the bubble closes in time t. This quantity 
can be resolved in experiments by measuring fiuorescence 
correlations of a tagged bubble Si. On plug ging Eq. 
into Eq. dM]), we get 



C(xo,t) = l- 



r(Ci + l/2,x;^/20 

r(Ci + 1/2) ' 



(25) 



where r(s,x) = t^~^ exp{—t)dt is the upper incom- 
plete gamma function. This result agrees with that re- 
ported in [1^. It is easy to derive the following asymp- 
totic behaviors of C(xo,t): In the limit t — > 00, one has 

M 



C{xo,t) 



(4) 



2-\Ci-|-l/2 



-Ci-1/2 



(Ci + l/2)r(Ci + 1/2) 
while, in the limit < ^ 0, one has 

C„2/9^Cl-l/2 

Cfxn « 1 - ^ ' ^1/2~C, -xl/2t 

^ °' ^ r(Ci + l/2) ^ 



(26) 



(27) 



B. Distribution of the area till the first-passage 
time: P{A\xo) 

The area A = J^^ x{t')dt' under the random motion 
([T7| can tell us about the readiness of the bubble to react. 
Here, the motion starts at xq and continues in time till 
the first-passage time. Note that the quantity A is not 
a geometric area, rather its units are length x time. To 
compute the related pdf, we substitute U{xo) = xq in 
Eq. (Ull), 



1 d^Qjxo) 

2 dxl 



Ci dQixo) 



Xo 



dxo 



-pxoQixo) = 0, 



which has the general solution [3J|, 

Q(xo) = x[;^+'/' \AiMiz) + A2J-,{iz) 



(28) 



(29) 



Here, z ^ (2/3)^/2^ Xq^^ u = (2Ci + l)/3, J^(x) is 
the Bessel function of the first kind, and Ai and A2 are 
arbitrary constants. Using J^{x) = i'^ Iv{—ix), where 
Ijjjc) is the modified Bessel function of the first kind 
[3l|, gives 



Q(xo) - x'^'+^f^ \BiI,{z) + S2/_.(^) 



(30) 



where Bi and B2 are arbitrary constants. Since for large 
X, /±i/(x) ~ /\/2ttx to satisfy the boundary con- 
dition Q{xo ^ 00) = 0, we must have Bi = — i?2- 
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To satisiy Q{xo 0) = 1, wc note that as xq —J- 0, 
one has I^{x) « {x/2Y/r{iy + 1) [HI, which yields 
B2 = r(l - i/)(v^/3)'^. We thus get 



Qixo) 



(I)' 



r(i 



r(^) 



I-,iz)~I,iz) 



(31) 



where, in obtaining the last relation, we have used the 
identities, Ki,{x) ~ (7r/2)[J^(x) — Iu{x)]/ sin(i'7r) and 



31| . On taking inverse Laplace 



exp(-2xg/9A) 



T{v)T(l~v) = IT / sin(z/7r) 

transform of ([5T|) , we obtain the desired pdf [32[ 

2(2Ci + l)/32,2Ci + l 
P{A\X0) = 3(4Ci+2)/3r((2Ci + l)/3) yl(2Ci+4)/3 " 

(32) 

This expression nicely reproduces the numerical results 
obtained by simulating the Langcvin equation ([3|) under 
the full potential with xq < Xch and C2 < 0, see the mid- 
dle panel of Fig. [3l The inset demonstrates an expected 
disagreement for larger bubbles with xq ^ Xch- 



C. Distribution of the maximum before the 
first-passage time: P{M) 

How large can the bubble grow before it closes, assum- 
ing an initial opening of xq? This question is of interest 
in estimating the efhciency of processes that can occur in- 
side big loops only. The relevant measure is provided by 
P{M), the pdf of the maximum bubble size M before its 
first closure, given that xq € [0, M]. We obtain this pdf 
by following closely the procedure of [H, [1^ ■ We first 
define a related function. Let q^xq) be the probability 
that the motion described by Eq. p7|) exits the interval 
[0, M] for the first time through the origin. Thus, q{xo) is 
the cumulative probability that the maximum before the 
first-passage time is < M. It is evident that this func- 
tion satisfies two boundary conditions: (i) q(0) = 1, and 
(ii) q{M) — 0. Denoting by 0Ar(Aa;) the distribution 
function of a small displacement Ax in time Ar — > 0, 
it follows from the Markovian property of the dynamics 
p7)) that 



q{xo + Aa;)0Ar(Aa;)d(Ax). 



(33) 



On Taylor expanding q{xo + Ax) and averaging over 
Ax = -iCi/xo)AT + mAT, using (^(0)) = 0, (^^(O)) = 
1/At, we get, to leading order in Ar, the equation 
[il/2)id^qixo)/dxl) - {Ci/x„)idqixo)/dxo)]AT = 0. For 
arbitrary Ar, wc obtain 

1 d'^q{xo) Ci dq{xo) _ ^ ^^^^ 



2 dxl 



0. 



xo dxo 

Solving this equation with the above mentioned bound- 
ary conditions, we get 

qixo)^!-^^'^' . (35) 



The pdf of interest is obtained by differentiating q{xo) 
with respect to M, 



P{M) 



(2Ci 



l)xf^^\ 



M > Xq. 



(36) 



In Fig. [3] (bottom) , we compare this result with numeri- 
cal simulations using the full potential. We observe good 
agreement for xq < Xch]C2 < 0, and an expected dis- 
agreement for Xq ^ Xch (inset). 



D. Joint pdf of the maximum M and the 
corresponding time tm before the first-passage time: 

P{M, tm) 

To compute the joint probability distribution of the 
maximum bubble size M and the time tm at which 
the maximum occurs before closure, we adopt the PD 
method, see Section Hi Dl and [23| . We split a path evolv- 
ing under (|17p into two parts, before and after the time 
t,n, with respective weights Wl and Wr. Due to the 
Markovian property, the weight W of the whole path is 
given by the product of the weights for the two split parts, 
W = Wl X Wr. The weight Wr can be obtained from 
Eq. ([55]) . Recall that Wr is the weight of a path that 
starts at M— e at time t = tm and exits the interval [0, M] 
for the first time through the origin. On the other hand, 
q{x{)) in Eq. p5p is the probability for a path starting at 
Xo e [0, Af] to exit the interval for the first time through 
the origin. We thus deduce that Wr, = q{M — e), or. 



Wr^\ 



(M-e) 



_ ,^(2Cl+l) 



(2Ci 



l)e 



M 



+ 0(e2). (37) 



The second equality is derived by assuming e to be in- 
finitesimal. The weight Wl is obtained from Eq. (fT5| by 
substituting C2 = 0, 



Wl (X 



M 



J2 e-^-*" 



i^piM - e)^Ppixo). (38) 



Here, ipp and Ep are the solutions of the eigenequation, 

1 d^ cl 



2dx2 2x2 



■0 = Elf)] < X < M, 



(39) 



subject to the condition tpi^x 
solution of this equation is 



i>p{x) = AVxJq {^y 



0, M) = 0. The general 



(40) 



where Ja(x) and Ya{x) are the Bcssel functions of or- 
der a o f the fir st and second kind, respectively, and 
a = iymcf [131 ■ Note that Ci is real, thus a > 0. 
Since for x ^ 0, Y^{x) « -(r(a)/7r)(2/x)" [U, we de- 
mand that _B = for satisfying ■i/'(x = 0) = 0. The 
other boundary condition results in the discrete eigen- 
values Ep such that 



where Unn denotes 
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FIG. 3: (Color online) Numerical simulations below the melt- 
ing temperature, T < Tm, for the pdf of the first-passage 
time (top) , the pdf of the area till the first-passage time (cen- 
ter) and the pdf of the maximum size till the first-passage 
time (bottom). The parameters are Ci = 1 and C2 = (red 
square), Ci = f and C2 = —0.1 (blue dots). The initial bub- 
ble size is 2:0 = 2 for all the cases. The analytic results in 
the small bubble approximation, (|22p . (|32}, (|36p appear in 
black continuous lines. The three insets compare the analytic 
small bubble results (black continuous line) with numerical 
simulations using C\ = 1 and C2 = —0.5 (purple +). 



the p-th zero of Ja{x). The constant A is determined by 
requiring ippi^) to be normalized. On using the identity, 

CdppJa{UapP/a)Ja{UaqP/a) = (5p,g(a^/2) [J^+i (Wap)]^ 

pl| . we finally get 



'2x 



M\Ja + l{Uap)\ 

and the probability 

\gi 2^(M-e)xo 
M-J A'P 

Uap{M 



J a 



M 



(41) 



Wt. oc 



^ [Ja+l{UapW 



M 



J. 



■A 



M J 



(42) 



Next, we evaluate Wl to leading order in e by Tay- 
lor expanding J a (^Vi^AM-J^ and also using the result 

J'aiuap) = —Ja+i{uap), which foUows from the following 
identity: J'^{z) = ^Ja{z) - Ja+i{z) We finally get 



Wl cx e 



2x. 



Ci+1/2 00 



p=\ 



-u^pt„/(2M^) 



fUapXo\ 

""^"V M ) 



+ 0(e2). 



(43) 



The probability P{M,tm; of the whole path is the 
product of ([43|) and ([38|) with a normalization con- 
stant i?(e), which is determined by requiring that 
lim^^o P{M,t,n;e)dt„, -> P{M), where P{M) is 
given by Eq. (pG]). 



(44) 



PiM,t,„;e) = Bie)WLWR. 



Using the identity J2^=i JaiuapX^/M)/ {uapJa+i{uap)) = 
x^/(2Af") for < xo < M [3^, we get B{e) = 

^ (1^)^-^^-'/', where C = 2C, + I - 
Substituting for B{e) in Eq. and taking the limit 

e — > 0, we get the desired probability. 



P(M, t™) = (2Ci + l) 

CXD 



= 1 



2^ 



(45) 



For the free Brownian motion, taking Ci — 0, this result 
reduces to that derived in |23||. 



IV. LARGE BUBBLE DYNAMICS 

We study here the dynamics of large bubbles of size 
X > Xc\i, see Eq. ([S]). In this limit, one can neglect the 
term —Ci/x in the Langevin equation ([3]), and study the 
dynamics dictated by 



dx 
'dt 



c2+m. 



(46) 
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This equation describes a one-dimensional random walk, 
x(t), in the presence of a constant drift, C2. The prob- 
ability distribution to find a bubble of size x at time t 
is a shifted Gaussian P{x,t) ~ g-(a:-xo-C2t)"/2t^ where 
xo is the initial size of the bubble. Below the melting 
temperature (C2 < 0), the bubble shrinks and eventually 
disappears in time. Above T^, when C2 > 0, the bub- 
ble grows in size as time passes. For large C2 > 0, the 
first-passage time diverges. Thus, our analysis below of 
various first-passage functionals is valid for large bubbles 
at temperatures below as well as above T„i, in the 
latter case while conditioned on a finite tf. 

For obtaining the pdfs, P{tf\xo) and P{A\xo), we 
adopt the BFP method. The relevant differential equa- 
tion is obtained from Eq. (|9]) by substituting Ci = 0, 



1 d^Qjxo) 

2 dx^ 



dxQ 



pUixo)Q{xo)^0, (47) 



with the boundary conditions (i) Q{xn — > 00) = 0, and 
(ii) Q{xo — >■ 0) = 1. We also derive analytical results for 
P{M) and P{M,tm)- 



A. First-passage time distribution: P{tf\xo) 

This distribution has already been investigated in (sBj 
in a different context. We thus omit the details of the 
calculation, but include the results for the sake of com- 
pleteness of our presentation. The procedure involves 
solving Eq. (|47)) with U{xo) = f under the boundary 
conditions, then taking inverse Laplace transform of the 
solution to yield 



Pitflxo) 



1 Xq 



cxp 



ixQ+C2tf)^ 

2tf 



(48) 



The moments (tj) may be obtained by using the identity, 
x'^-ie-^/^-'^^da; = 2{l3/jY/'^K^{2y/^) for Re(/3) > 
and Re(7) > 0, where KJx) is the modified Bessel 
function of the second kind [301 ■ One gets 



(4) = ^^(|^)'(xo|C2|)i/2e-'^-»X,_V2(|t^2|xo). 



(49) 

e ■", the mean first- 
-2C2X0 



/C2 for C2 > 0, 



Noting that Ki/2{x) = ^7r/(2x) 
passage time is given by (tf) — xqc 
and by (tf) = X0/IC2I for C2 < 0. The survival proba- 
bility, defined in Eq. (fM)) . is given by 

C{xo,t) = 1 - 4Le-^^^^« /*i7^/2e-^o/(2t/)-c,=*//2rfi^, 
V27r Jo 

(50) 

For C2 < (T < Tm) and for large t, one has 
Cixo,t) « f - (C2Xo)/(2V^)e-c^^-«--o/(2t)[r(-f/2) - 

r(-l/2,C|t/2)j. Using the resuh that r(s,x) 




0.8 
0.7 

0.6 h 

^ 0.5 

o 

^ 0.4 - 



0.3 
0.2 
0.1 




tf 



• ■ • 1 



• 

■ 




X 


X • ■ 


• ■ 




■ 

X • 

■ 




•m 








2 


3 4 5 




A 




0.01 



FIG. 4: (Color online) Results from numerical simulations be- 
low the melting temperature for the pdf of the first-passage 
time (top), the pdf of the area till the first-passage time (cen- 
ter), and the pdf of the maximum size till the first-passage 
time pdf (bottom). The parameters are Ci = (red square), 
Ci — 1 (blue dots), and Ci = 4 (purple cross). The initial 
bubble size is a;o = 5 and C2 ~ —4 in all the cases. The 
analytic results in the large bubble approximation appear in 
black continuous lines. In the middle panel, we do not make 
a comparison with the analytic function for P(yl|a::o) since its 
explicit form is known only in the limit of small and large 
values of A. 
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FIG. 5: (Color online) Results from numerical simulations 
below and above Tm, showing the first-passage time pdf (top), 
the pdf of the area till the first-passage time (center), and the 
pdf of the maximum till the first-passage time (bottom). The 
parameters are C2 = —2 (red filled square), C2 = — 1 (blue 
filled circle), C2 = —0.5 (purple cross), C2 = 0.1 (light blue 
filled triangle), and C2 = 0.5 (black circle). The initial bubble 
size is a;o = 2 and Ci = 1 in all cases. 



B. Distribution of the area till the first-passage 
time: P{A\xo) 

This quantity can be obtained by solving Eq. (j47p with 
U{xo) = xo with the appropriate boundary conditions, 
then deri ving the inverse Laplace transform of the solu- 
tion, see [35l. Is?! for details. In particular, one obtains 
the following two limiting behaviors of the distribution 
P{A\xo): For A — > c>o, one has 



PiA\xo) 



e-^^-«sinh(|C2|a;o) {2^^^ f\C2\Y^^ 



A 



1/2 



In the opposite A ^ Umit, one gets j37 | 



P{A\xo) 



32/3r(l/3) A^/3 



(51) 



(52) 



Note the distinct asymptotic forms in Eqs. (ISTI) and (|52p . 
While the latter demonstrates a behavior similar to that 
observed in the small bubble dynamics [cf. Eq. ([5^ ]. the 
former predicts a different scaling behavior. 

Results from numerical simulations are displayed in 
Figs. |3]-[5] For a fixed value of C2 < 0, on increasing 
Ci, Fig. 13] shows the narrowing of P{A\xo) and the dis- 
placement of its center towards smaller A values, thereby 
reflecting the increased importance of the entropy term 
in the free energy function. On fixing Ci and on increas- 
ing C2 from negative to (small) positive values, the area 
pdf develops an increasing contribution at large A val- 
ues, thereby hinting at the onset of large bubbles and 
long first-passage times when T > Tm- 



00 



[31[, we get C{xo,t) 

y27^(xo/C|)t-3/2exp [-(xo - |C2|t)V(20]- 

Figure m compares the analytic result for P{tf\xo) with 
numerical simulations obtained by considering the full 
potential. Since, for xq > Xch with C2 > 0, many tra- 
jectories have diverging tf,wc performed here numerical 
simulations only below the melting temperature, so that 
C2 < 0. In this case, with increasing IC2I, the effective 
drift velocity towards bubble closure increases and bub- 
bles quickly disappear in time. The top panel demon- 
strates that the analytical prediction (|48)) agrees with 
simulation results for a;o > Xch]C2 < 0. Figure [5] further 
includes results from numerical simulations by consider- 
ing the full bubble potential, displaying the behaviors 
both below and above the melting temperature. Upon 
increasing C2 from negative values (i.e., for T < Tm) to 
positive values (i.e., for T > T^), one notes that (i) the 
center of the pdf P{tf\xo) is displaced to longer times, 
and that (ii) the bubble lifetime is significantly enhanced 
(see Fig. [51 top panel) . This corroborates with the phys- 
ical picture that with increasing temperature, a bubble 
takes a longer time to disappear. 



C. Distribution of the maximum before the 
first-passage time: P{M) 

The procedure here proceeds as in Section fill CI One 
first finds the cumulative probability q{xo), which satis- 
fies m, m 



1 d'^q{xo) dq{xo) _ 

' L^2 — : — U. 



2 dx'^ 



dxn 



(53) 



With the boundary conditions (i) q{Q) = 1 and (ii) 
q{M) = 0, this gives the solution 



qixo) 



smh[C2(M - Xq)] 
sinh(C2M) 



(54) 



The desired pdf is obtained as the derivative of q{xo) 
with respect to M, 



P{M) 



C2e-g^"° sinh(C2Xo) 
sinh2(C2M) 



A/>xo. (55) 



11 



This result has been obtained in [23|,|35| in a different con- 
text. Figure 0] compares this form with numerical results 
at temperatures below T,„ (C2 < 0). We observe a good 
agreement for initial bubble sizes satisfying xq > Xch- 
Figure [5] further displays numerical results, both below 
Tm and above it, by adopting the full bubble potential. 



D. Joint pdf of the maximum M and the 
corresponding time tm before the first-passage time: 

PiM,tm) 

This distribution has been studied in [23l | for the ran- 
dom motion (|46|) with C2 < by employing the PD 
method described in Section III Dl To compute P{M,tm) 
for general C2 , we follow the discussion of Section IIII Dl 
and split a representative path into two parts, before 
and after tm- We compute the weight of each part sep- 
arately. The weight Wr of the path after t„i is given by 
Wr = q{M - e). Using Eq. ([M]), and then taking e to 
be infinitesimally small, we get 



,-C2{M-e) 



sinh(C2e) 



G 



26 



-C2M 



sinh(C2A/) 



'smh{C2M) 



(56) 

The weight Wl of the part before is obtained from 
Eq. (Uni) by substituting Ci = 0, 



Wl oc eC2(*^-^-"«) ^ e-^''*"^p(Af - e)Vp(2;o). (57) 

p=i 

Here, ipp{x) and Ep{x) are the eigenfunctions and ener- 
gies of the Hamiltonian (|T3l) with the potential. 



Vix) 



^ if < x < M, 
00 if X = 0, a; = M. 



(58) 



The norm alized eigenfunctions are easily obtained as 
ippix) = y^2/M sm{pTTX / M) with the corresponding en- 
ergies Ep = C|/2-hp27r2/(2M2). Substituting this in Eq. 
([571). we get 



Wl oc ee 



C2(A/-a:o)-C|t„./2 



27r 



00 

5](-l)^+Vsin(^)e-P^-^*-/(2^^^)+0(e2). 



(59) 



The last equation, together with Eq. (j55|) . gives the total 
probability. 



P{M,tm;e) = B{e)WLWR. 



(60) 



The normalization constant B{e) is determined by re- 
quiring that lim^^o P{M,tm',i)dtm — > P{M), given 
in Eq. 1^. With X;^i(-1)^~ Vsin(px)/(p2 + a^) = 
(7r/2) sinh(a2;)/ sinh(a7r) for — tt < a < tt [30I |. one gets 
B{e) = l/2e^. Using this expression for B{e) in Eq. (pn|) 
and taking the limit e ^ 0, we get the desired pdf, 



P(M,t„ 



sinh(C2A/) M2 



X ^(-l)f+Vsin 
p=i 



\ M ) ' 



(61) 



The marginal distribution P{t„i) can be obtained by in- 
tegrating P{M,tjn) over M from xq to infinity. The 
large-im and small- tm asymptotic behaviors of P{tm) for 
C2 < are discussed in [2J. 
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Quantities 


Results for small bubble 


Results for large bubble 


Pitflxo) 


^f]^'-^'^ ■ tf ^OO 


^ ^y^/^^-{xa+C2tfy /2tf 

^ ij3/2g_C.Xoe-C^//2 . tf^ oo 

^ g-(xo+C.t,)V2t/ . t^. _^ 


C(xo,<) 


~ r(Ci+i/2) e ' , t ^ u 


[(xo-ICalt)^ 1 


PiA\xo) 


^ ^-i(Ci+2)g-2xi5/9A 
_^-|(Ci+2). 
_g-2.3/M. ^^Q 


_ ^-4/3g-2.g/9A . ^ ^ 


P{M) 


- Af-2Ci-2 


1 

^ sinh^(C2M) 




OO ,,2 J. /Q™r2 

,,,r_(c+3) e "q,p'™'''^" j- /M„,pi:o\ 
p=l 

where C = 2Ci + 1 - ^1 + 4C2/2 


p=l 



Table 1: Scaling behavior of the probability distribution functions of various Brownian functionals calculated in this 
work for small and large DNA bubbles. 



CONCLUSIONS 



appropriate boundary conditions. 



In this paper, we derived probability distribution func- 
tions of various Brownian functionals associated with a 
random walk model for DNA bubble dynamics at tem- 
peratures below, at, and above the denaturation tem- 
perature. Based on the backward Fokker-Planck method 
discussed in [2^ . we derived (i) the first-passage time 
distribution P{tf\xo), providing information about the 
bubble lifetime, (ii) the distribution P{A\xo), of the area 
A covered by the random walk till the first-passage time, 
measuring the bubble reactivity to processes within the 
DNA, and (iii) the distribution P{M), of the maximum 
bubble size M before bubble closure, all conditioned on 
an initial bubble of size xq {xq £ [0,M]). (iv) The joint 
probability distribution P{M,tm) of the maximum bub- 
ble size M and the time i,„ of its occurrence before the 
first passage time was also obtained by employing the 
Feynman-Kac path integral formulation. The advantage 
of the elegant methods adopted here is that they produce 
results on various functionals by making proper choices 
of a single term in a parent differential equation with 



We considered separately the dynamics of small and 
large bubbles. Analytical results for the pdfs at each limit 
nicely agree with Langevin simulations. Our analysis re- 
veals different nontrivial scaling behaviors of P(ty|a;o), 
P(^|a;o), P{M) and P{M,tm), as summarized in Ta- 
ble 1. The scaling exponents are characterized either 
by the entropic parameter Ci, or by the base-pair dis- 
sociation parameter C2. These quantities may thus be 
estimated experimentally by using fluorescence correla- 
tion spectroscopy Q to measure, e.g., the maximum size 
distribution P{M) for small and large bubbles separately. 



We expect our results to be useful in quantifying chem- 
ical processes within DNA, for example, protein binding 
to single-stranded DNA, and for developing a deeper un- 
derstanding of polymer dynamics. It is of interest to ex- 
tend our study and consider loop-loop interactions [27[, 
and the effects of disorder and heterogeneity in predict- 
ing the kinetics of specific processes within DNA bubbles 



13 



Acknowledgments 

DS and MB acknowledge support from the Connaught 
fund and from NSERC. SG thanks Amir Bar, Raphael 
Chctrite, Ori Hirschbcrg, Satya N. Majumdar and David 



Mukamel for fruitful discussions and suggestions, and 
gratefully acknowledges the Israel Science Foundation 
(ISF) for supporting his research at the Weizmann In- 
stitute of Science. 



J. D. Watson and F. H. C. Crick, Cold Spring Harb. 
Symp. Quant. Biol. 18, 123 (1953). [22 
A. Kornberg and T. A. Baker, DNA Replication Second [23 
Edition (University Science Books, Sausalito, California, 
2005). [24 
M. D. Frank-Kamenetskii, Phys. Rep. 288, 13 (1997). [25 
R. M. Wartell and A. S. Benight, Phys. Rep. 126, 67 
(1985). [26 
D. Poland and H. A. Scheraga, Theory of Helix- Coil 
Transitions in Bio-Polymers (Academic Press, New [27 
York, 1970). 

C. Danilowicz, Y. Kafri, R. S. Conroy, V. W. Coljee, [28 
J. Weeks, and M. Prentiss, Phys. Rev. Lett. 93, 078101 
(2004). 

O. Krichevsky and G. Bonnet, Rep. Prog. Phys. 65, 251 
(2002). 

G. Altan-Bonnet, A. Libchaber, and O. Krichevsky, [29 
Phys. Rev. Lett. 90, 138101 (2003). 

J. Phys.: Condens. Matter, Special Section on DNA melt- [30 
ing 21 (2009). 

R. Metzler, T. Ambjornsson, A. Hanke, and H. C. [31 
Fogedby, J. Phys.: Condens. Matter 21, 034111 (2009), 
and references therein. 

H. Kunz, R. Livi, and A. Siito, J. Stat. Mech.: Theory [32 
Exp. P06004 (2007). 

T. Ambjornsson, S. K. Banik, M. A. Lomholt, and R. 

Metzler, Phys. Rev. E 75, 021908 (2007). 

S. K. Banik, T. Ambjornsson, and R. Metzler, Europhys. 

Lett. 71, 852 (2005). [33 

A. Hanke, and R. Metzler, J. Phys. A: Math. Gen. 36, [34 

L473 (2003). 

A. Bar, Y. Kafri, and D. Mukamel, Phys. Rev. Lett. 98, [35 
038103 (2007). 

A. Bar, Y. Kafri, and D. Mukamel, J. Phys.: Condens. [36 
Matter 21, 034110 (2009). 

T. Dauxois, M. Peyrard, and A. R. Bishop, Phys. Rev. 
E 47, 684 (1993). [37 
S. Srivastava and Y. Singh, Europhys. Lett. 85, 38001 
(2009). [38 
H. C. Fogedby and R. Metzler, Phys. Rev. Lett. 98, 
070601 (2007); Phys. Rev. E 76, 061915 (2007). [39 
L.-A. Wu, S. S. Wu, and D. Segal, Phys. Rev. E 79, 
061901 (2009). [40 
Note that we refer to the functionals here as "Brownian 
functionals" , although the dynamics involves both the [41 
random delta-correlated force and a deterministic com- 
ponent corresponding to the bubble potential, see Eq. 



©■ 

S. N. Majumdar, Current Science 89, 2076 (2005). 

J. Randon- Furling and S. N. Majumdar, J. Stat. Mech.: 

Theory Exp. P10008 (2007). 

M. Kac, Trans. Am. Math. Soc. 65, 1 (1949). 

S. N. Majumdar and M. J. Kearney, Phys. Rev. E 76, 

031130 (2007). 

P. L. Krapivsky, S. N. Majumdar and A. Rosso, J. Phys. 
A: Math. Theor. 43, 315001 (2010). 

Y. Kafri, D. Mukamel, and L. Politi, Physica A 306, 39 
(2002). 

As discussed in [S^l , one is required to impose an infinites- 
imal cut-off e, or else an infinite number of crossings of 
the level M (immediately after the first contact at time 
tm) would make it impossible to ensure that the path 
stays below Af at times subsequent to tm- 
A. C. Brahka and D. M. Heyes, Phys. Rev. E 58, 2611 
(1998). 

I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, 
Series, and Products (Academic Press, New York, 1965). 
G. B. Arfken and H. J. Weber, Mathematical Methods 
for Physicists Fourth Edition (Academic Press, London, 
1995). 

Tables of Integral Transforms, Volume I, Based, in part, 
on notes left by Harry Bateman and compiled by the staff 
of the Bateman Manuscript Project, edited by A. Erdelyi, 
M. F. Oberhettinger and F. G. Tricomi (McGraw-HiU, 
New York, 1954). 

A. J. Bray, Phys. Rev. E 62, 103 (2000). 

Introduction to Bessel Functions, F. Bowman (Dover, 

New York, 1958). 

M. J. Kearney and S. N. Majumdar, J. Phys. A: Math. 
Gen. 38, 4097 (2005). 

A. P. Prudnikov, Yu. A. Brychkov and O. L Marichev In- 
tegrals and Series: Volume 2 Special Functions (Gordon 
and Breach Science Publishers, Netherlands, 1998). 

M. J. Kearney, S. N. Majumdar, and R. J. Martin, J. 
Phys. A: Math. Theor. 40, F863 (2007). 

B. Coluzzi and E. Yeramian, Eur. Phys. J. B 56, 349 
(2007). 

T. Ambjornsson, S. K. Banik, O. Krichevsky, and R. 

Metzler, Phys. Rev. Lett. 97, 128105 (2006). 

T. Ambjornsson, S. K. Banik, O. Krichevsky, and R. 

Metzler, Biophys. J. 92, 2674 (2007). 

M. Zoh, Phys. Rev. E 81, 051910 (2010). 



